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Abstract. - We examine the temperature dependence of thermal conductivity of one dimensional 
nonlinear (anharmonic) lattices with and without on-site potential. It is found from computer 
simulation that the heat conductivity depends on temperature via the strength of nonlinearity. 
Based on this correlation, we make a conjecture in the effective phonon theory that the mean- 
free-path of the effective phonon is inversely proportional to the strength of nonlinearity. We 
demonstrate analytically and numerically that the temperature behavior of the heat conductivity 
k oc 1/T is not universal for ID harmonic lattices with a small nonlinear perturbation. The 
computer simulations of temperature dependence of heat conductivity in general ID nonlinear 
lattices are in good agreements with our theoretic predictions. Possible experimental test is 
discussed. 
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The role of nonlinearity in the dynamics of one dimen- 
sional (ID) nonlinear (anharmonic) lattices has attracted 
attention for many decades, for instance, the ergodicity 
problem in the Fermi-Pasta-Ulam (FPU) chains inducted 
by nonlinear interactions [1,2], solitons in nonlinear partial 
differential equations [3] and discrete breathers in nonlin- 
ear lattices [4,5], to name just a few. Recent years have 
witnessed increasing studies on the role of nonlinearity in 
heat conduction in low dimensional systems (See Ref. [6] 
and the references therein). The fundamental question 
is whether the non-integrability or chaos is an essential 
or a sufficient condition for the heat conduction to obey 
the Fourier law. From computer simulations, it is found 
that in some ID nonlinear lattices such as the Frenkel- 
Kontorova (FK) model and the <p A model, the heat conduc- 
tivity is size independent [7-9] . Thus the heat conduction 
in these models obeys the Fourier's law. This transport 
behavior is called normal heat conduction. Whereas in 
some other nonlinear lattices such as the FPU and alike 
models, the heat conduction exhibits anomalous behav- 
ior [10], namely the heat conductivity k diverges with the 
system size N as K ~ N s . Studies in past years have fo- 
cused on the physical origin and the value of the divergent 
exponent S [10-14]. It is found that the anomalous heat 
conduction is due to the anomalous diffusion and a quanti- 



tative connection between them has been established [15]. 
Most recently, a very general effective phonon theory has 
been proposed to describe the normal and anomalous heat 
conduction under the same framework [16]. 

As a by-product of the study of heat conduction in 
low dimensional systems, nonlinearity (anharmonicity) 
has been found very useful in controlling heat flow. Be- 
cause of the nonlinearity, the lattice vibrational spectrum 
depends on temperature. This property has been used to 
design the thermal rectifiers/diodes [17-22] and thermal 
transistors [23]. Inspired by the two segment theoretical 
models [18,19], Chang et al [24] has built the first solid 
state thermal rectifier with a single walled carbon nan- 
otube and boron nitride nanotube, which indicates the 
opening of a new research field - controlling heat flow at 
microscopic level by using nonlinearity. 

However, in contrast to the size-dependence, the tem- 
perature dependence of heat conductivity and its relation- 
ship with the nonlinearity have not yet been studied sys- 
tematically even though this problem is very fundamental 
and very important from the experimental point of view. 
In fact, for low dimensional nanoscale systems such as nan- 
otube and nanowires etc, to measure the temperature de- 
pendence of thermal conductivity for a fixed length sample 
is much easier than that for the size dependence of ther- 
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mal conductivity for a fixed temperature. The difficulty 
for measuring size-dependent thermal conductivities lies 
in the fabrication of low dimensional systems. One can- 
not guarantee that the systems (nanotube/nanowires) of 
different system sizes fabricated separately are identical. 
However, performing the experiments at different temper- 
atures for a fixed length is feasible. 

In this Letter, we study the temperature dependence of 
heat conductivity and its connection with nonlinearity in 
several representative ID nonlinear lattices with and with- 
out on-site potential, namely, the FK model, the </> 4 model, 
and the FPU models. We will show that the temperature 
dependence of heat conductivity cannot be understood by 
the phenomenological phonon collision theory [25] in har- 
monic lattice with a small nonlinear perturbation. We 
find via computer simulations that the heat conductivity 
is actually inversely proportional to the strength of non- 
linearity. More interestingly, if we make a conjecture in 
the effective phonon theory that the mean-free-path of ef- 
fective phonon is inversely proportional to the strength of 
nonlinearity, we can give a rather satisfactory explanation 
for the temperature dependence of heat conductivities at 
both weak and strong coupling regimes in ID nonlinear 
lattices such as the FK and the FPU-/? models consis- 
tently. We argue that the results presented in the current 
paper can be tested by nano scale experiments. 

Without loss of generality and simplicity, we consider a 
ID nonlinear lattice whose Hamiltonian reads, 



N 



^-+V(x i ,x i+1 ) + U(x i ) 



(1) 



where Xi is the displacement from the equilibrium position 
of the ith particle. In the following study we consider both 
cases with on-site potential, such as the FK model and the 
</> 4 model, and the cases without on-site potential such as 
the FPU-a, the FPU-/? model and the FPU-a/? model. 

The Frenkel-Kontorova model In the FK model V = 
\{xi — a^+i) 2 , and U = , 2 ^ s (1 — cos27nEi). In this Let- 
ter, we take parameter K = 2tt. The FK model is very 
special since it reduces to a harmonic lattice with a small 
nonlinear perturbation at both low and high temperature 
regimes. According to the phonon collision theory, the col- 
lision frequency of a given phonon should be proportional 
to the number of phonons with which it can collide. Since 
the phonon distribution function is proportional to T in 
classical regime, the total number of excited phonons is 
proportional to T. Hence the mean-free-path of phonon 
I oc 1/T. At this regime, the heat capacity and velocity 
of phonons are all constants. From the Debye formula, 
one gets the heat conductivity, k oc I. Therefore the heat 
conductivity, n oc 1/T universally for the harmonic lattice 
with a small nonlinear perturbation. 

We compute the heat conductivity of the FK model as 
a function of temperature with fixed boundary conditions 
and the Nose-Hoover heat baths. The results are shown 
in Fig.l in double logarithmic scales. At the high temper- 



ature regime, the on-site potential is considered as a small 
perturbation, thus the FK model is suitable for the in- 
terpretation of phonon collision theory. However the heat 
conductivity in this regime is not inversely proportional 
to temperature n oc 1/T as the phonon collision theory 
predicts, it increases almost linearly with the tempera- 
ture k oc T! This striking contradictory suggests that the 
phonon collision theory may not be suitable to describe 
the classical heat transport. 

To explain the temperature dependence of heat conduc- 
tivity in ID nonlinear lattices, we introduce a new quan- 
tity to describe the strength of nonlinearity We define a 
dimensionless nonlinearity e as a ratio between the aver- 
age of nonlinear potential energy and the total potential 
energy which consists of linear and nonlinear potential en- 
ergy: 

, \{En) 



{Ei + E n ) 



< e < 1, 



(2) 



where (•) is an ensemble average. For the harmonic lattice 
with a small nonlinear perturbation, the average of nonlin- 
ear potential energy (E n ) is negligible compared to the av- 
erage of linear potential energy (Ei), thus the dimension- 
less nonlinearity can be approximated by e s=s ^4 ^ , where 
the average of the linear potential energy (Ei) = NksT/2 
at the harmonic limit according to the theory of energy 
cquipartition [26]. 

At the high temperature regime, the average of nonlin- 
ear on-site potential energy (E n ) is negligible compared 
to the average of harmonic inter particle potential energy 
(Ei). The dimensionless nonlinearity can be expressed as 



( Ei£ (1-008 27^ )) 



N 



, (3) 



where the average of the on-site potential energy of 
any particle, (1 — cos 2itXi), has been assumed to be a 
constant in the interval (0, 2) at the high temperature 
regime. The average of inter particle potential energy 

(EiLi 3 ( x i - x i+if) = Nk B T/2. Therefore, the dimen- 
sionless nonlinearity at high temperature regime has the 
asymptotic behavior: e oc 1/T. Compared with the com- 
puter simulation it is not difficult to find the correlation 
between the heat conductivity and the dimensionless non- 
linearity at the high temperature regime for the FK model: 



k oc - oc T. 



(4) 



This correlation strongly suggests that the mean-free-path 
of phonons should be proportional to the inverse of the 
dimensionless nonlinearity, I oc 1/e, rather than to the 
inverse of temperature, / oc 1/T, as the phonon collision 
theory predicts. 

As for the nonlinear lattice, the notion of phonon is not 
valid any more, one should invoke the effective phonon 
theory [16]. Recently, we find out that the normal and 
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anomalous heat conduction can be treated with effective 
phonon theory [16]. For a general ID nonlinear lattice 
with 



V( Xi ,x i+1 ) = ^> X : +l) , U(xi) = J>A (5) 

s=2 



the spectrum of effective phonons is u)\ = a (1^+7), where 
lo = 2y / 52sinfc/2 is the phonon dispersion of harmonic 
lattice, and the temperature dependent coefficients a and 
7 are defined as 



a = 



(Lite - x i+ i) 2 ) 



7 



(6) 



The heat conductivity can be expressed by the modified 
Debye formula: 

c f 27T 

k=— P(k)v(k)l(k)dk, (7) 

where c is the specific heat, P(k) the normalized power 
spectrum of total heat flux, v(k) — duj^/dk the velocity 
of the effective phonon, and l(k) the mean- free-path of 
effective phonon. l(k) = v(k)r(k), where r(fc) is the effec- 
tive phonon relaxation time which is proportional to the 
quasi-period of effective phonons, r(k) = A27r/a)fe, and the 
dimensionless A only depends on temperature. 

Here we conjecture that the mean-free-path of the effec- 
tive phonon is inversely proportional to the dimensionless 
nonlinearity, 

Acxi. (8) 

In the harmonic limit e — > 0, the phonons should have 
infinite mean-free-path which is exactly implied by this 
conjecture. The temperature dependence of heat conduc- 
tivity can be expressed as 



-^—V oc -l—V, 

2tt e 



(9) 



where 



Jo 



P ( k )-T^2 — ^TTT? cos2 77 dfc - 
; K+7) 3/2 2 



The low and high temperature regimes are of special in- 
terest, in which the heat conductivity exhibits asymptotic 
behaviors. The transport properties of general nonlinear 
lattices are mostly determined by the leading terms of 
potential energy perturbed by another potential energy 
with the first order approximation. For example, the di- 
mensionless FPU-/3 lattice has quadratic and quartic inter 
particle potential energy 



i=l 



Pi 1, 



Xi+1) 2 + \{xi - x i+1 ) 4 



• (10) 



At low temperature regime, the leading term of poten- 
tial energy is the quadratic, and the quartic term is re- 
garded as a small perturbation. Whereas at high temper- 
ature regime, the leading term is the quartic one, and the 



quadratic term is a small perturbation. At these temper- 
ature regimes we can assume that P(k) is temperature in- 
dependent since the topology of the Hamiltonian doesn't 
change with temperature. Furthermore, the unchanged 
topology of the Hamiltonian also implies that the specific 
heat c is almost temperature independent. In the classical 
case, c = (H)/NT where (H) is the ensemble average of 
total energy at temperature T. The specific heat c = 1 at 
classic harmonic lattice if we take the Boltzman constant 
fc_B to be unity. We take the FPU-/? lattice as example. 
At low temperature regime the FPU-/3 lattice reduces to 
a harmonic lattice, so c « 1. At the high temperature 
regime, the leading order of potential energy is the quartic 
potential energy. Due to the general equipartition theory, 

NT = (£i Siff ) » (£ite - Xi+i) 4 )- Then the specific 
heat cw (^PV^) /NT + - x l+i r) /NT = 

3/4 where we have ignored the harmonic potential at high 
temperature regime. The specific heat c is just a constant 
in both temperature regions. This is true for general ID 
nonlinear lattices, thus we can write, 



K (X 



V. 



(11) 



Therefore, the heat conductivity depends on temperature 
via the dimensionless nonlinearity e. The two coefficients 
a and 7 can be calculated from statistical mechanics. 

For the FK model, the coefficient a = 1 exactly. At the 
high temperature regime, the dimensionless nonlinearity 
e cx 1/T as we have discussed above, thus the heat con- 
ductivity depends on temperature as 



k oc TV. 



(12) 



Since 7 goes to as the temperature goes to infinity, 7 
should decrease with the increase of temperature. As a 
result, V should increase with temperature T. Therefore, 
the heat conductivity should increase faster than a linear 
dependence of T. This is what we observe in numerical 
results shown in Fig.l, where n oc T 130 at the high tem- 
perature regime. 

At the low temperature regime, the on-site potential can 
be expanded in a Taylor series: 



U( Xi ) 



1 

2^ 



{2-KXi) 2 (27TX;) 4 (2-rrXi) 



2! 



4! 



+ 



6! 



(13) 



The linear potential energy consists of the harmonic in- 
ter particle potential energy and the harmonic on-site po- 

+ 2tt|| . The 
Nk B T/2. The 



tential energy, E t = £\ =1 \{ Xi - x i+1 ) 

equipartition of energy still holds, (Ei) = 
first order nonlinearity is estimated by the average of 
quartic on-site potential over the linear potential energy, 

6 ~ (J2iLi -h (27 4? ) )/(^) oc T- However, the particles 
are not stable in the Hamiltonian with this first order 
nonlinear approximation because of the minus sign be- 
fore the quartic on-site potential. At this low tempera- 
ture regime, the system is most likely to approach a sta- 
ble state with the second order nonlinear approximation. 
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Fig. 1: (Color online)Heat conductivity k vs temperature, T, 
for the FK model with lattice length N = 500. The best fit 
line is k oc j 1-1 - 9 ™- 1 a t the low temperature regime and is 
k oc J 11 - 30 * - 03 a t the high temperature regime. The dotted line 
in the middle is the turning point estimated from the Hamil- 
tonian, To = 2K / (2-ir) 2 , which is 1/-7T, in the case of K = 2ir in 
this paper. 



The proper estimation of the dimensionless nonlinearity 

should be proportional to (J2iLi 2IF ^§T~)/(^i) ^ ''■ 
Therefore, the heat conductivity should depend on tem- 
perature as k oc 1 /T 2 if we neglect the slight temperature 
dependence caused by the coefficient 7. This is also in 
a good agreement with the numerical result k oc 1/T 1,9 
in Fig.l. The slight difference might come from the first 
order unstable nonlinear potential and the temperature 
dependence caused by the coefficient 7 which has been 
neglected in this situation. 

From Fig.f it is clear that the heat conductivity de- 
creases at low temperature region and increases at high 
temperature region. The turning point can be roughly 
estimated from the Hamiltonian itself. Here we use the 
expression of the specific heat again 



c = 



K 



-l) 2 + — 



7 (f — cos 27^ 



NT 



(14) 

By definition, (J2iPl/ 2 ) = NT/2. At this turning point, 
we expect the ensemble average of inter particle potential 
and the on-site potential to be the same value 



K 



V( 2 *) 



■(1- cos 27^)). (15) 



The temperature of the turning point thus can be ex- 
pressed as 



T = 



2K 



(E^-co^tt^)) 2K 



1/2 



N 



(16) 



where we have used the approximations c rs 1 and 
(2,(1 — cos 2-7^)) w N/2. With K = 2ir in our numeri- 
cal calculation, we have Tq — 1/tt, which is shown as the 



dotted line in the middle of Fig.l. It is in a good agree- 
ment with the numerical calculations. 

The cf> 4 model Another commonly studied model is the 
(f) 4, model which has an on-site potential U(x) — \x A . The 
coefficient a = 1. At the low temperature regime, the <f> 
model is topologically similar to the FK model. The non- 
linear quartic on-site potential in the (j) 4 model guarantees 
the particles in stable states, thus the heat conductivity 
depends on temperature as n oc V/T. V decreases with 
temperature because 7 increases with temperature. Thus 
the heat conductivity should decrease faster than k oc 1 /T 
which is consistent with the numerical results in Ref. [9] 
where k oc 1/T 135 . The high temperature regime of the 
<j) model is not in our consideration because it has no 
energy transport at high temperature limit which is also 
called anti-continuous or anti-integrable limit [4]. 

The FPU-P model In the FPU-/3 model, U(x) = 0, and 



Xi+i) + j(xi — x i+ i) . We take k = 1 and 
(3=1. For the general ID nonlinear lattices without on- 
site potential such as the FPU-/3 model, they have total 
different dynamical properties at low temperature (weak 
coupling) regime and high temperature (strong coupling) 
regime. The leading order potential energy is quadratic 
potential at low temperature regime and quartic potential 
at high temperature regime. The coefficient 7 = for gen- 
eral ID nonlinear lattices without on-site potential. This 
property allows us to predict the dependence of heat con- 
ductivity more precisely since V is a constant. In this 
case, 

(17) 

At low temperature regime, the coefficient a ~ 1 by 
definition and the heat conductivity is just proportional to 
1/e. The dimensionless nonlinearity for the FPU-/3 model 
is e w |(xi - Xi+i) 4 )/ (J2i ~ Xi+i) 2 ) oc T. Thus 
the heat conductivity k oc 1/T which is consistent with 
the numerical results in Ref. [27] and in the Fig. 2. 

At high temperature regime for the FPU-/? model, the 
linear potential energy is negligible compared to the non- 
linear quartic potential energy. The dimensionless non- 



linearity is e = 



1. The 



(EiLi §(*i-xi + i) 2 +§(zi-x i+ i) 4 ) 

heat conductivity only depends on coefficient a at this 
region as 

k oc \fa. (18) 

The coefficient a can be exactly calculated from the for- 
mula in Ref. [16] 
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2 ^(4) 211/2 ji 



00. 



(19) 



Therefore, the heat conductivity at high temperature 
regime is n oc T 1 / 4 . This temperature behavior is consis- 
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Fig. 2: (Color-online) (a) Heat conductivity kvsT for the FPU 
lattices with TV = 2000. The solid circles are numerical results 
of the FPU-/9 model. Open circles are that of the FPU-a model 
with absolute cubic potential where data have been shifted ver- 
tically for comparison. The inset is high temperature behavior 
of the FPU-a model which is actually calculated with param- 
eters k = and a* = 1. k = is equivalent to the high 
temperature limit. The lower panel is the enlargement of low 
temperature part. The additional semi circles are from the 
FPU-q/3 model with cubic and quartic interparticle potential. 



tent with previous results in Ref. [27] and our numerical 
results in Fig. 2. 

The FPU-a model We consider the FPU-a model in 
which V = \{xi — x i+ i) 2 + 2L- \xi — x i+ i\ 3 . We take k = 1 
and a* = 1. Here the absolute cubic inter particle po- 
tential is used to prevent the particles from escaping to 
infinity which occurs in the usual FPU-a model. At the 
low temperature regime of the FK model, we have met 
the case that the heat conductivity doesn't decrease with 
temperature as k oc 1 /T which is believed to be universal 
in solid state physics. Here we present another case that 
the heat conductivity is not proportional to 1/T at the 
low temperature regime. 

Fig. 2 shows the numerical results of heat conduc- 
tivity of this modified FPU-a model as a function 
of temperature. At the low temperature regime, 



The dimensionless nonlincarity e rs 
iX /(£i5(zi-*i+i) 2 ) « VT. Thus 



oc 1/VT. 

3 \x% 

the prediction derived from our conjecture k oc 1/e oc 
1 / \T at low temperature regime is in good agreement 
with the numerical results. 

At high temperature regime, the coefficient a is 



/ °Ve V 2 + 3 J /T d0 



Jo°V 



Jo°V 



T 1/3 . (20) 



The heat conductivity should depend on temperature as 
k oc y/a oc T 1 / 6 . The numerical results in Fig. 2 con- 
firm this prediction. For the general realistic solids, the 
inter particle potential can always be expanded as a poly- 
nomial form. However, the first order nonlinear poten- 
tial is cubic where particles will escape to infinity. Parti- 
cles are more likely to approach the stable states governed 
by the next order nonlinear potential: quartic potential. 
This might be the reason that the realistic solids have 
the universal behavior (k oc experimentally at this 

regime. To check this in ID we calculate the temperature 
dependence of heat conductivity of the FPU-a/3 model 
which has both cubic and quartic inter particle poten- 
tial besides the harmonic inter particle potential, namely, 



V{xi,x t+ i) = 1 {x l -x i+ i)- + \{x i 



) 3 + f(* 



i) 4 



(with fc = l, a* = 1, (3 = 1) and find the temperature de- 
pendence of k is more like that of FPU-/? model (lower 
panel of Fig.2). 

Finally we discuss the possible experiment to verify our 
conjecture. From the computer simulations we can see 
that eventually the nonlincarity will cause the heat con- 
ductivity to bend upward as the temperature increases. 
It is worthwhile to point out that the <?i 4 model does not 
show this property because it cannot model realistic solids 
in the high temperature regime. It is not physical to make 
the substrate on-site potential larger than the interparti- 
cle potential which binds the atoms together. The reason 
we cannot observe this kind of temperature dependence of 
heat conductivity in realistic solids at room temperature 
is that the strength of nonlincarity in them is too weak. 
The realistic solids just start melting before the nonlinear- 
ity dominates the transport. However, with nano technol- 
ogy we can build the ID sample attached to a substrate 
to demonstrate this counter-intuitive temperature behav- 
ior of heat conductivity at room temperature caused by 
nonlinearity. Now we consider the proposed lattice model 
with real physical units: 



N r 



p z 1 

— + -mu 2 (x i+1 - x{f + U{Xi) 
2m 2 



(21) 



where m is the mass of atom, ui is the oscillation fre- 
quency, and U(xi) the weak substrate on-site potential 
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with strength A which can be controlled by experiment. 
As we have discussed above, this on-site potential U{x{) 
does not need to be the FK like potential. U(xt) can 
be any kind of on-site potential with strong nonlinearity. 
Here the "strong" means the linear part inside U(xi) is 
small compared with the nonlinear part. If the magni- 
tude of this on-site coupling strength is A, then using the 
same analysis we have used for the FK model, the heat 
conductivity should increase with temperature after 



To = - — (Kelvin), 
kb 



(22) 



and we predict the increase is larger than the linear in- 
crease with temperature. To observe this classical behav- 
ior explicitly, we must lie outside of the quantum regime 
where the increase of heat conductivity is mostly caused 
by the increase of specific heat with temperature. Finally, 
we must emphasize that the key point for this experiment 
is the substrate on-site potential with "strong" nonlinear- 
ity. 

In summary, we have found via computer simulation the 
strong correlation between the temperature dependence of 
heat conductivity and the strength of nonlinearity in ID 
nonlinear lattices. In the effective phonon theory of heat 
conduction, we conjecture that the mean-free-path of ef- 
fective phonon is inversely proportional to the strength of 
dimcnsionless nonlinearity. The predictions from this con- 
jecture arc consistent with the numerical results in several 
very general ID nonlinear lattices with and without on- 
site potential. We found the temperature dependence of 
heat conductivity k oc 1/T is not universal for the model 
of harmonic lattice with a small nonlinear perturbation. 
The temperature dependence of heat conductivities in the 
commonly studied FK and FPU models are explained by 
our theory at both weak and strong coupling regimes con- 
sistently. We have also proposed possible experiment to 
verify our conjecture. 
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